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Abstract 

^r*j ■ Formation and evolution of general relativistic collapse inhomogeneity is investigated. 

It is shown that a rarefaction wave forms at initial collapse stages and propagates inside 
i \ from the cloud boundary to its center. The focusing time of the rarefaction wave is found. 

■ In massive clouds, the rarefaction wave focusing time equals the free-fall time. The collapse 

of such clouds leads to black hole formation. In low-mass clouds, the rarefaction wave 
ly-j | focusing time is less than the free-fall time. After focusing, the collapse of such clouds 

' becomes fully inhomogeneous and can be sufficiently decelerated by the pressure gradient. 

o 

1 Introduction 

The study of general relativistic gravitational collapse is a complex problem because of math- 
ematical difficulties appearing when solving the corresponding non-linear partial differential 
equations even in spherical symmetry pQ . Only a few exact solutions are known for any special 
cases More realistic models for gravitational collapse can be investigated using numerical 
• i—i ! simulations (see, e.g., [3]). 

An initially homologous collapsing cloud with zero pressure (dust cloud) remains homologous 
during the collapse. But contraction of an initially homologous cloud with non-zero pressure 
develops with time in a sufficiently non-homologous way because a pressure gradient must 
necessarily develop at the cloud surface {3J. Even if the cloud is initially surrounded by an 
ambient gas of the same pressure as the cloud, the further contraction breaks this equilibrium 
[Hj. A rarefaction wave propagates in from the cloud surface to its center with sound speed. 
The rarefaction front divides the cloud into two parts (01 QUI- The pressure gradient becomes 
zero in the internal region. However, some inhomogeneity forms in the external region. This 
dynamical inhomogeneity can evolve in a self-similar manner at last stages of contraction [SI |§] . 

The dynamics of rarefaction wave has been studied for non-relativistic (protostellar) collapse 
of an isothermal cloud [ID] . The effects of rotation ^T] and magnetic field [7] have been taken 
into account. In this paper, a general- relativistic analogue of this problem is considered under 
the condition of spherical symmetry. 

In the next section, the basic equations are described. In Sec. 3, the flow structure in the 
inner homogeneously collapsing core is considered. The equation of rarefaction front propaga- 
tion is solved in Sec. 4. In Sec. 5, the radial null geodesies are investigated. Finally, Sec. 6 
presents the main results and conclusions. 
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2 Basic equations 



Let us work in a system of spherical coordinates (t, R, 6, if) moving at each point with the gas 
located at that point (comoving or Lagrangian coordinates). The interior space-time geometry 
for a collapsing cloud in spherical symmetry can be described by the metric of a "collapsing 
star" (see, e.g., [T2]): 

ds 2 = e 2 *c 2 dt 2 -e 2A dR 2 -r 2 dn 2 , (1) 

where dVt 2 = d6 2 + sin 2 Odip 2 . Here $, A, r are some functions of the radial coordinate R and 
the time t to be determined by the Einstein equations. Note that the form of the line element 
(JIJ is invariant under the transformation: R! = R'(R), t' = t'(t). 

In the comoving reference frame, the components of the four- velocity of gas are u° = e - *, 
u 1 = u 2 = u 3 = 0. Therefore, for a perfect fluid, the stress-energy tensor has following non-zero 
components: Tq = e, T\ = T 2 — T| = — P, where e is the gas energy density (gas energy per 
unit proper rest volume) and P is the gas pressure. 

Using these relations, we can obtain the following Einstein field equations: 

(2) 



(3) 




f' : r' , ,r 



= A--$'-, (4) 

tyt iy tyt 

where the prime means a derivative with respect to the radial coordinate R and the overdot 
means a derivative with respect to the time t. The equations of gas motion give 

e f P' 

-^ = -a-2-, = -<*>'■ (5) 

e + P r e + P 

The system (J2JEJ) is still not complete. We must add any equation of state to close this 
system. In this paper, we assume the following equation of state: 

P=(K-l)e, (6) 

where k = const (1 < k < 2). From this equation it follows that the sound speed in the gas is 
a = \J n — lc. 

Eq. © can correspond to two reasonable astrophysical situations 8 . In the first case 
(protocluster and protogalactic clouds) this equation can be considered as the equation of state 
of a non-relativistic isothermal ideal gas when adiabatic heating and cooling is neglected and 
therefore the internal energy is negligible compared to the rest energy. In the second case 
(presupernova cores), Eq. ((Hj) can be considered as an extreme relativistic limit of an adiabatic 
(or polytropic) equation of state when the rest energy density is negligible relative to the internal 
energy. 



3 Homogeneously collapsing core 

Inside the internal homogeneous core of a collapsing cloud, the pressure gradient is zero. From 
the second equation in (J3J) it follows that $' = and therefore $ = $(£)• Thus we can define 
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another time variable T = T(t) such that $ = 0. It is clear that the time T coincides with the 
proper time of the gas. In addition, suppose that the radial coordinate R satisfies the relation 
r(T,R) = R(l — T](T)), where rj(T) is a function of the proper time T such that < 77 < 1. 
The value 77 = 1 corresponds to free-fall time (without the pressure gradient) for the considered 
homogeneous region. 

From Eq. (J1J we get the following solution: 



1—7] 



A = In 

y/l + fW 

where f(R) is some function. Combining (J7J) and (|5)i. we find 

e 



(7) 



3k ' 



(1-77) 

where eo is initial energy density. 

Now we can use Eqs. (J7J), (jHJ) to transform Eq. (J3J) to the following form: 



8ttG 



rj 



f 



r 



C 2 (l-77) 2 R 2 (l-T]) 2 R(l-T)) 2 



(9) 



Since the left-hand side of this equation is a function of the proper time T, we see that / = foR 2 , 
where f is some constant. Taking into account the initial values 77(0) = 1 and 77(0) = 0, we 
can find: f = — 87rGe /(3c 4 ). 

The invariant line element JIJ can be written as 



/ 



ds 2 = c 2 dT 2 



dR 2 



\ 



1 



87rGeo 
3c 2 



+ R 2 dtt 2 



R 2 



(10) 



Note that this metric corresponds to the well-known Friedmann solution for a closed homoge- 
neous isotropic universe with uniform pressure. 

Now we can combine the expression for energy density (JHJ with to obtain an equation 
for the function 77 (T): 

l87cGe r 1 1/2 



drj 
dT 



3c 2 



1 



(1 - ?7) 3k 



3k-2 



1 
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This equation can be solved analytically in the special case k = 4/3. For an arbitrary parameter 
ac, Eq. (|11|) can be solved numerically. 



4 Propagation of a rarefaction wave 

The rarefaction wave front propagates in collapsing gas with a sound speed. The location 
R r f(T) of the rarefaction wave front satisfies the equation: 
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Figure 1: Functions Rrfiv) i n the case K — 5/3 for various values of parameter 



a. 



After some manipulations, we can transform Eq. (jl2j) to the form: 

dr rf V«- 1(1 - rf)^ y/l - a 2 r 2 f 



d V ay/l - (1 

where r r f = R r f/Ro, Ro is the initial cloud radius, 

[r, 



v) 



3k-2 



(13) 



a 




, Rn 





rin 



2GM _ 8nGe R 3 Q 
c 2 ~ 3d ' 



(14) 



The parameter a controls the general relativity effect in the problem. In the case of very 
small a, the space-time curvature is negligible, and the dynamics of the rarefaction wave can 
be considered as non-relativistic. In the opposite limit, as a — > 1, the problem is highly 
general- relativistic. 

The solution of Eq. ()13)) can be written as: 



1 • ( ■ l^f^l M . 

r r f = — sin I arcsm a H — — arccos(l — rj) 



a 



2 



3k - 2 



(15) 



The solutions obtained, for the case k = 5/3 and various a, are shown in Fig. 1. 

The rarefaction wave front reaches the center of the cloud at a focusing time T*. The value 
of T* can be determined by an equivalent value = ^(T*). There are two possible regimes for 
propagation of rarefaction waves in general relativistic collapsing clouds. 

For sufficiently large a > a*, where 



7Ti//t — 1 

a* = sin- — . 

3n — 2 



(16) 



we have 77* = 1. This result means that the rarefaction wave front reaches the cloud center at a 
moment of big crunch when the central homogeneous core collapses to space-time singularity. 
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For relatively small values of parameter a < a* the rarefaction wave front reaches the cloud 
center before the formation of a central singularity. The focusing time is determined by the 
value 

-2- ( 3k ~ 2 \ 

77* = 1 — cos 3 *- 2 — . arcsina . (17) 
\2v« - 1 / 

In such a cloud, after focusing of the rarefaction wave, a strong pressure gradient can be formed. 
Hence the further collapse can be sufficiently decelerated and can even pass to slower quasi- 
hydrostatic contraction with formation of some stable object at final stages of cloud evolution 
(e.g., a neutron star). Note that, for a fixed initial radius Rq, a small value of the parameter a 
corresponds to a small initial mass of the cloud. The simplest expression for the focusing time 
can be obtained in case k = 4/3: 



T* = J — — — sin (V3 arcsin a 
V 87rGe v 

The critical value a* is a function of the parameter k. As k —>■ we have the asymptotic 
relation a* = 7Ta/ k — 1. In particular, in the case of a non-relativistic isothermal gas, the 
obtained criterion coincides with the non-relativistic one (see [TTHllll IT]). The maximum value 
of a* is reached at K — 4/3. At this point, a* = sin7r^/3/6 ~ 0.79, and therefore the initial 
cloud radius is R ~ 1.61/L. 



5 Null radial geodesies 

We can probe the space-time geometry in the internal homogeneous region by emitting outgoing 
light rays. The equation for outgoing radial null geodesies in the inner homogeneous region is 

After some simple manipulations we can transform this equation to the following form: 




dR / 3c 4 v /; V 3c 4 

1 (20) 



drj V 87rGe — (1 — ^) 3 «- 2 

For null geodesies that are outgoing from the cloud center, we have the following solution: 

Rq flf , 3k-2 3k-2 -| 'l 

R = — sin< arcsin(l — r^o) 2 — arcsin(l — m 2 >. (21) 

a 1 3/t — 2 L ' J J 



Here 7] is an initial value of the parameter r] for a given outgoing light ray. 

Fig. 2 shows the space-time geometry in a collapsing cloud with the parameters k = 5/3 
and a = 0.85 in terms of the coordinates 77 and R. The heavy solid line corresponds to the 
rarefaction wave front. The arrowed lines denote outgoing radial light rays. Each arrow points 
the direction of light ray propagation. Some outgoing rays can reach the rarefaction wave front 
and then get out to the outer inhomogeneous region of the collapsing cloud. Other rays cannot 
reach the rarefaction wave front because they fall back to the cloud center at 77 = 1. Every light 
ray emitted at any point of the shaded area cannot escape the internal homogeneous region. 
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Figure 2: Space-time geometry of a collapsing cloud in case k = 5/3 and a = 0.85. Lines with 
arrows denote outgoing radial light rays with various values of the initial parameter r] . 



In other words, we can conclude that an event horizon forms in the considered cloud at final 
stages of its collapse. 

Therefore, for clouds with a < a*, there are no effects that can retard the collapse, and 
some part of matter necessarily falls under an event horizon with the formation of a black hole. 
The obtained criterion for black hole formation in collapsing clouds can be reformulated in 
another form: Mq > M*, where the critical mass M* = 0.5alc 2 Ro/G depends on the initial 
cloud radius Rq. Note that in the 1 the value of the critical mass M* can be very small 

because in such a cloud the critical parameter a* is very close to zero. On the other hand, an 
event horizon can appear in a collapsing cloud after the rarefaction wave focusing. Hence the 
real critical mass can even be less than M*. 

6 Conclusion 

The basic results of this work are: 

1. The problem of inhomogeneity of general relativistic collapse of an initially homogeneous 
cloud surrounded by an ambient gas with the same pressure is investigated. It is shown that, 
at initial collapse stages, a rarefaction wave forms at the cloud boundary and then propagates 
to the cloud center. The front of this rarefaction wave divides the cloud into two parts. In the 
inner homogeneous region, the space-time metric is like the Friedmann one. In the outer region 
behind the rarefaction front, a strong inhomogeneity forms. 

2. The rarefaction wave front propagates in the collapsing gas with the sound speed. The 
characteristic equation of the rarefaction front propagation has been solved. The focusing time 
of the rarefaction wave has been found. 
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3. An analysis reveals two possible regimes of rarefaction wave dynamics. In very massive 
clouds, the rarefaction wave focusing occurs at the central singularity formation time (big 
crunch). The collapse of such clouds leads inevitably to black hole formation. In a cloud with 
a relatively small mass, rarefaction wave focusing takes place before the big crunch. After the 
rarefaction wave focusing, the collapse of such a cloud can pass to quasi-static contraction with 
further formation of some stable object. 
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